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Abstract: We identify various universal contributions to the entanglement entropy 
for massive free fields. As well as the 'area' terms found in [1], we find other geometric 
contributions of the form discussed in [2]. We also compute analogous contributions 
for a strongly coupled field theory using the AdS/CFT correspondence. In this case, 
we find the results for strong and weak coupling do not agree. 
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1 Introduction 

Entanglement entropy has emerged as a topic of interest in a wide variety range of 
research areas ranging from condensed matter physics [3] to quantum gravity [4, 5]. 
In the context of quantum field theory (QFT), when one considers the entanglement 
between two regions, 1 one finds that the entanglement entropy is UV divergent because 
of short range correlations in the vicinity of the 'entangling surface' E separating the 
two regions. If the calculation is regulated with a short distance cut-off 5, the leading 
contributions for a QFT in D spacetime dimensions generically take the form 

where each of the coefficients C2k involves an integration over the boundary S. For 
example, the leading term then yields the famous 'area law' result with c 2 oc [4]. 
Unfortunately the coefficients appearing in these power law divergent terms above are 
sensitive to the details of the UV regulator — for further discussion, see section 4. 



1 These are spatial regions on a fixed Cauchy surface. 
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However, certain subleading contributions can reveal universal data describing the 
character and/or the state of the underlying QFT. A well-known example of such a 
contribution arises for conformal field theories (CFT's) in an even number of spacetime 
dimensions [6-12]. Here, in calculating the entanglement entropy, one typically finds 
a logarthmic contribution \og(5) where the coefficient is some linear combination of 
the central charges appearing in the trace anomaly of the CFT. The precise linear 
combination is again determined by an integral of various geometric factors over the 
entangling surface. For a four- dimensional CFT, this universal contribution takes the 
form [9] 



In 



C abcd h ac h bd - Kl b Kl a + l -Kl a K\ b 



(1.2) 

where a and c are the usual central charges which appear in the trace anomaly [13]. 
The various geometric factors include: h ao , 71%, K l ah , respectively, the induced metric, 
intrinsic scalar curvature and extrinsic curvature of E; C a b c d, the (pull-back of the) 
Weyl curvature of the background geometry; and L, some characteristic scale in the 
geometry. 

A similar class of universal contributions were identified in [1] for massive QFT's. 2 
In particular, considering free massive scalar fields, the following universal contribution 
to the entanglement entropy was found 



I 7r> As vnP 2 log(m<f) for even D , 
univ ~ 1 lD ATm D - 2 for odd D . 



where m is the mass of the scalar and Ay, is the area of the entangling surface. For the 
free scalar theory, the numerical coefficient 7 D is given by [1] 



{-) D > 2 

D-2 



6(47r)^-r(D/2) 

'Yd, scalar \ d — i V / 



12(4tt)- 2 2- 2 T(D/2) 



for even D 
for odd D . 



Motivated by these free field results, ref. [2] began a study of analogous terms 
for strongly coupled field theories using holographic techniques [8, 16]. This approach 
allows one to study the effect of perturbing a UV fixed point QFT by introducing 
general relevant operators, beyond simple mass terms. The results reveal a broad class 
of new universal contributions to entanglement entropy, which can be schematically 



^Related results were found previously in [14] and [15]. 
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represented as 



S^ = ,(D,n) [d°->*Vh [Mature 



mP 2 2n log mS for even .D , 

m D-2-2n for odd ^ 



(1.5) 

where m is the mass scale appearing in the coupling of the relevant operator. The 
schematic geometric factor denotes a combination of both the background and extrinsic 
curvatures with a combined dimension 2n. Note that < n < (D — 2)/2 in these 
expressions. Hence in general, the entanglement entropy contains a family of universal 
contributions that includes the 'area' terms in eq. (1.3) at n = and purely geometric 
terms, analogous to those in eq. (1.2), at n = (D — 2)/2. Ref. [2] focussed on universal 
contributions proportional to \og(8), however, note that in general these terms (1.5) can 
appear for either even or odd D if the operator dimension of the relevant deformation is 
chosen appropriately. 3 The analysis of [2] also readily extends to show the appearance 
of new universal terms generalizing the cut-off independent terms in eq. (1.3) for odd D. 
An important distinction is, however, that generally these cut-off independent terms 
depend on the underlying state of the boundary theory while the logarithmic terms are 
state independent. 

One of our objectives in the following is to extend the analysis of [1] to reveal the 
new curvature terms appearing in eq. (1.5). This is easily accomplished by performing 
the free field calculations on a curved background. In particular, ref. [1] work with 
a 'waveguide' geometry, M. 2 x I D ~ 2 , where I is a finite interval with either Dirichlet 
or Neumann boundary conditions imposed at the endpoints. In section 2, we perform 
analogous calculations on a 'spherical waveguide', M 2 x § D ~ 2 . Our analysis also provides 
a number of other interesting extensions of that in [1]. As well as considering the 
entanglement entropy, we also calculate the Renyi entropy, which is another useful 
measure of entanglement [17, 18]. In section 2.1, we consider a free massive scalar field 
but we also include the possibility of the curvature coupling |£i?0 2 . In section 2.2, we 
extend the analysis to consider a free massive fermion. For the fermions, we find that 
the coefficients of the universal 'area' terms in eq. (1.3) become 



(-) D 2 /2 

6(271-)^ T(D/2) 



for even D 



Id, fermion = ^L2J ^D, scalar = \ p-i V / 

for odd D , 



(-)' 

U=T 



I2yj2{2-K)— T- T(D/2) 



3 For example, choosing an operator with A = (D + 2)/2, the corresponding coupling would take 
the form gm^ D ~ 2 ^ 2 where g is a dimensionless coefficient. In this case, an 'area' term with n = 
appears for all D > 3 with 7 oc g 2 . 
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as was noted previously in [19]. Above, /2\ denotes the integer part of D/2. Our 
analysis in these sections also allows us to identify a particular curvature contribution 
(1.5) to the entanglement entropy as 

f ,n 9 rr m/ ,s f rn D ~^ log(mS) for even D > 4 , ,„ 

s„„, v = 7D y^vAK( ft )x| mD _ 4 foroddD >v M 

where IZ(h) is the Ricci scalar of the metric induced on the entangling surface. For a 
free scalar with a curvature coupling, the new numerical coefficient jr> is given by 

Id, scalar = ~ ( £ — lD,scalar (1-8) 



2 V 6, 

where ^o^caiar is precisely the coefficient appearing in eq. (1.4). For a massive free 
fermion, this coefficient can be expressed as 

D-2 ^ 

T-D, fermion ~2A 'YD, fermion ( 1-9) 

where j D ,fermion is given by eq. (1.6). 

In section 3, we turn to a holographic calculation of entanglement entropy for the 
strongly coupled M = 2* gauge theory, a massive deformation of the celebrated M = 4 
super- Yang-Mills theory [20, 21]. While similar calculations already appear in [2], the 
details of the boundary mass terms and their translation to the dual gravity theory 
are precisely understood in this well-studied framework [22-24]. Hence, we are able to 
compare the 'area' contribution (1.3) at strong coupling from holography to the weak 
coupling results, which combine eqs. (1.4) and (1.6). Our final conclusion is that the 
strong coupling result does not match the corresponding contribution (1.3) found at 
weak coupling! 

We conclude the paper with a brief discussion of our results and future directions 
in section 4. This is followed by two appendices presenting various technical results: 
Appendix A explicitly demonstrates the validity of the separation of variables (2.11) 
used in section 2.2 for spin-| fields. Finally, Appendix B repeats the calculations for 
section 2 for a 'hyperbolic' waveguide IR 2 x HP 0-2 . 



2 Renyi entropy on a spherical waveguide 

In general to define the entanglement or Renyi entropies of some quantum system, we 
begin by dividing the degrees of freedom into two subsets, A and A, and calculate 
the reduced density matrix pa = Tr^ (p) where p describes the global state of the 
system. In a QFT, this is typically realized by beginning with a Cauchy surface in 
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a fixed background and introducing of an 'entangling surface' £, which divides this 
surface in two separate regions, again, denoted A and A. The reduced density matrix 
Pa is then given by integrating over all field configurations in the region A. Given this 
density matrix, the entanglement entropy is then defined by the standard von Neumann 
formula 

S EE = -Ti(p A \ogp A ) . (2.1) 
while the Renyi entropy is given by [17] 

S a = — ^- log Tr (p A a ) . (2.2) 
1 — a 

The latter is usually evaluated for (positive) integer values of a, in which case, eq. (2.2) 
involves a somewhat simpler calculation since it does not require evaluating the log- 
arithm of p A appearing in eq. (2.1). Further, if the result of S a can be continued to 
real values of a, (as will be possible in the following,) the entanglement entropy can be 
recovered as the limit: S EE = lim^i S a . 

The trace required in eq. (2.2) has a standard path integral representation, e.g., 
[25, 26], 

Tr(p A a ) = Z a /(Z 1 ) a , (2.3) 

where a is again assumed to take integer values for the moment. Implicitly, the first 
step here is to Wick rotate the background geometry to Euclidean time t E = it. Then 
Z a is the partition function of the QFT evaluated on an a-fold cover of the Euclidean 
background where a cut is introduced throughout region A on the Cauchy surface, 
which we denote t E = 0. At the cut, the fields on the fc'th sheet are joined to the fields 
on the (fc+l)'th sheet when approaching from t E — > 0" and to those on the (&-l)'th 
sheet when approaching from t E — > + . Hence Z a is the partition function evaluated 
on a singular covering geometry with an angular excess of 2n(a — 1) at the entangling 
surface E. The factors of the standard partition function Z\ appear in the denominator 
of eq. (2.3) to ensure that the density matrix is properly normalized with Tr (p A ) = 1. 
Given these partition functions, the Renyi entropy (2.2) becomes 

S a = kg^-alogZi {2A) 
1 — a 

Before proceeding with our explicit calculations, let us introduce the following 
shorthand notation for simplicity 

d = D-2. (2.5) 

Now following [1], our calculations here will focus on free fields in a waveguide 
geometry. In particular, throughout this section, the background geometry will take 
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the form M. = IR 2 x § d , where M. 2 is covered by Cartesian coordinates, t E and x. As 
in the above discussion, the Cauchy surface is selected by simply setting t E = in 
which case the resulting spatial slice has the waveguide geometry R x S d . This slice 
is then divided into two halves by choosing the entangling surface to be the sphere at 
x = 0. Implicitly we assume that the QFT on Ai is in its ground state and as described 
above, we are considering the reduced density matrix on the region A = {t E = 0, x > 0} 
resulting after integrating out the field degrees of freedom in A. 

To calculate the Renyi entropy as described above, we need to evaluate the partition 
function on Ai a , the a-fold cover of Ai. At this point, we note that the background 
geometry has a rotational symmetry in the plane around the point (t E , x) = (0, 0), which 
serves as our entangling surface. Hence in constructing Ai a , we are simply replacing the 
M. 2 component of Ai by a two-dimensional cone C a with an angular excess of 2ir(a — 1) 
at the origin, i.e., Ai a = C a x E> d . For later calculations, we explicitly write the metric 
on C a as 4 

ds 2 = a~ 2 (dr) 2 + r 2 (d9) 2 , (2.6) 

where r and 9 possessing the full radial and angular range < r < oo, < 9 < 2ir. 
Now this geometry has no distinguishing features which prefer integer values of a 
(apart from a = 1) and so from this point forward, we allow a to take any (positive) 
real value. That is, we are analytically continuing a already in the covering geometry 
[27] rather than first evaluating Z a for integer a and then analytically continuing. 
The rotational symmetry in the transverse space around the entangling surface is an 
essential ingredient for this geometric approach. 5 

One feature, which distinguishes our background here from that in [1], is that the 
cross-section of the waveguide § d is curved and hence our calculations of the entan- 
glement entropy below can reveal new universal contributions of the form given in 
eq. (1.5). After a few more preliminary remarks, we will consider a massive free scalar 
field in section 2.1, with the action 

5(0) = / dPxy/g\( ( V0) 2 + m 2 2 + ill 2 ) , (2.7) 

where TZ is the Ricci scalar of the background geometry. Hence we have included a non- 
minimal coupling to the curvature of the background here. Of course, with m = and 
£ = 4(£)~ 2 i) ) we have a conformal scalar field theory. In section 2.2, we also consider a 

4 After changing variables r — > r/a, 9 — > aO, one can readily conclude that the angular excess is 
given by 2ir(a — 1). However, we use the present coordinates (2.6) in the following, since the angular 
momentum operator takes the standard form in this representation of the cone. 

5 See [10] for further discussion. 
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massive free fermion field (with minimal coupling to the background geometry) , which 
becomes a conformal theory if m = 0. 

For either of the above classes of theories, the partition function is Gaussian and 
can be exactly evaluated using the heat kernel approach, e.g., [28] 



A2-KS POO 71 

p t 



where (t,x,y) with x, y G M. a is the heat kernel of the corresponding massless 
wave operator on A^ a . The trace of the heat kernel involves taking the limit of coin- 
cident points, i.e., y — > x, and integrating over the remaining position x. Of course, a 
trace is also taken over the spinor indices in the case of the spin-| field — see below. 
In the above expression and throughout the following, we use s = or | to indicate the 
scalar or fermion cases, respectively. Further, 6 is a short-distance scale introduced to 
regulate any potential UV divergences (as discussed in the introduction). Finally m s 
denotes the 'effective' mass of the field under study. For the fermion, we have simply 
m s= i = m, however, given the non-minimal coupling of the scalar in eq. (2.7), we have 



' 2 



-Lo = - 2 + ^^^, (2.9) 

where the second term comes from the curvature of the S d , i.e., 1Z{E> d ) = d(d - l)/R 2 
for a sphere of radius R. 

Notice that the curvature of the full background geometry contains a singularity 
at the tip of the cone, e.g., 

K(M a ) = TZ(E> d ) + K{C a ) = di * d ~ l) + 4tt(1 - a) S^(r) + ... (2.10) 

where r*G C a and the second term corresponds to the leading contribution from lZ(C a ) 
in an expansion near a ~ 1. To treat this singularity in a well defined way, we delete 
the point at the tip of C a and work on the space (C a — {0}) x S d . Of course, this 
means that appropriate boundary conditions must be imposed at the tip to make the 
wave operator self-adjoint. This then becomes the requirement that we must use only 
non-singular eigenf unctions when the heat kernel is constructed [29]. To simplify the 
notation we use the notation Ai a = C a x E> d to denote the punctured manifold in what 
follows. 

The wave operators are separable on the product manifold C a x S d and hence the 
heat kernel on A^„ can be expressed as the product of the two individual heat kernels 
on C a and E> d , i.e., 

K$ = (2-11) 
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where for brevity we have suppressed the arguments of the heat kernels here. Note 
that while this separation of variables is obvious in the case of the spin-0 field, it is 
less evident in the case of the spin-| field due to spinor structure of the heat kernel. 
Hence we show that separation of variables indeed holds in the latter case in Appendix 
A. Given eq. (2.11), one can write 

TrK$ a = r rrK% r ItK$, (2.12) 

where each trace on the right-hand side involves an integration over the corresponding 
component of the product manifold. In the case of spin-| field, there is also trace over 
spinor indices. Using the conventions established in Appendix A, we can regard the 
two traces on the right as also including a separate trace over the spinor spaces of the 
two component manifolds, C a and § d . 

In eq. (2.8), the possible UV divergences at t — > in the partition function were 
regulated by introducing the short- distance cut-off 5. However, we would also like to 
introduce a ^-function regularization [28] here since it readily allows us to identify the 
universal contributions to the Renyi entropy for general values of d. This approach will 
be applied in calculating to the Renyi entropy of the scalar field on C a x E> d in the next 
subsection. However, we also apply this regularization to produce general results for 
both scalars and fermions on the hyperbolic waveguide C a x M. d in Appendix B. 

In the (^-function approach, the partition function is regulated by shifting the power 
of t in eq. (2.8) 

i2irs roo r. 

logZ^ = e —5-^J o -^Trtf&e-**, (2 . 13) 

where now 5 appears to keep the whole expression dimensionless. Of course, after 
carrying out the integral over t, the regulator must be removed by taking the limit 
z —7- and suitably renormalizing the parameters of the theory to eliminate possible 
divergences in z. 

Now, as shown in [29, 30], the trace of the heat kernel on a cone depends on a only, 
therefore substituting eq. (2.12) into eq. (2.13) yields 

log Z« = — 5-^ TrJfW v{z) C^}{z) , (2.14) 
where the (^-function is defined as follows 

Expanding eq. (2.14) in the vicinity of z = then yields 



logZ^ = — TrK t 



2 



dz 



-4?(0)log5 2 + O(z) 

z=0 



(2.16) 
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where we rescaled 5 2 — > e _7 5 2 to absorb a term proportional to the Euler constant 
7. The pole term in the above expression must be removed by suitably renormalizing 
the field theory parameters. 6 The remaining contributions are precisely those which 
determine universal contributions to the Renyi entropy 



i2irs 



logZ« = — TrK t 



2 



dC 



dz 



C£ } (0)log^ 



(2.17) 



This expression is readily evaluated using results available in the literature [31, 32]. 
The ^-function on the hyperbolic space M d was computed for the spin-0 case in [31] 
and for spin-| case in [32]. For the scalars, the desired ^-functions on S d are then easily 
obtained from the hyperbolic ones using a formula given in [31] - see eq. (2.25). 

2.1 Renyi entropy for a massive scalar 

We begin here by evaluating the partition function Z a for the massive free scalar field 
theory described by the action (2.7). In this case, the heat kernel in eq. (2.8) corre- 
sponds to the D-dimensional scalar Laplacian on M. a . Separation of the variables leads 
to eq. (2.12) and further the heat kernel on C a is given by [33] 

Tri ^ (t) = 1^ (1 " ^ + a Tri ^ (t) • (2 ' 18) 
On the other hand, the scalar heat kernel on S d satisfies 

(-d t + AW) K$(t, x,y) = 0, x,yeS d (2.19) 

K$(0,x,y) = 5{x,y) , 

where is the scalar Laplacian on E> d . Of course, due to the rotational symmetry 
the heat kernel only depends on the arc-length between the two points on the sphere. 
Therefore for simplicity, we place one of the points at the north pole of the sphere. 
With this choice, the heat kernel becomes a function of only the azimuthal angle 9 and 
we may replace A^ by its radial part. 

To solve the resulting equation, we follow the prescription described in [32]. That is, 
we consider the intertwining operator, O = —(2ix sin 6)~ 1 dg, which relates the Laplacian 
on spheres of different dimensions, i.e., 

Afo = O^A { e d - 2) + , A { e d) = d 2 e + (d-l)cot8d e . (2.20) 



6 Appearance of the logarithmic term in cqs. (2.16) and (2.17) might seem misleading. Indeed, such 
terms are expected in the case of even d only. However, as follows from eq. (2.25) and eq. (B.4) for 
s = or eq. (B.25) for s = 1/2, Csd (0) = for odd d and so the log<5 term vanishes as expected. 
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Here R is the radius of both E> d and E> d ~ 2 . The overall constant factor in O is chosen 
such that O relates the delta functions on the two spheres. Hence, for even and odd 
dimensions, we have 

KSUt,6) = e^'* (-^—- e 9e) n K^(t,6) , (2.21) 
KjgUW) = e^l* (__L_Q,y K £\t,9) . (2.22) 

For d — 1, Kg\t) can be readily evaluated using the method of images. It is 
given by an infinite sum of the scalar heat kernels on R, which are shifted by integer 
multiples of 2ir with respect to each other to maintain periodic boundary conditions 
for the scalar field on a circle, namely 



oo 



Kg\t,B) = -== £ *— L *- L - ■ (2-23) 

V 47TT 



n=— oo 



For d = 2, the heat kernel can be constructed using spherical harmonics Yi m (9,<j)), 
which correspond to the orthonormal eigenf unctions of the Laplacian on a unit two- 
sphere. If one of the points is taken to the north pole, the result then simplifies to the 
following sum 

K * {t ' B) = 4^5> / + ^(cos^e- 1 ^ . (2.24) 

1=0 

We will also use the ^-function formulae given in eq. (2.17). As shown in [31], 
^-functions on §> d and M d for d > 3 are related and can be obtained from each other by 
means of complex contours. The final result in the case of spin-0 case reads 



$(z) = e m(2 - d/2) Vol(§ d ) 



iR 2 Z - de -i ^ d 2 r(c?/2) s . n{nz) (2 25) 



Vol(H d ) 7i d / 2 

f(ib - y)dy 

o [1 + (-iy e My-ib)] [( y _ ihf + b 2 Y 



with Re(z) < 1 



where denotes the scalar ^-function on H d , which is given by eqs. (B.4) and (B.6). 
Following the notation of [31], we also have the following definitions: 

b 2 = -tf m l + V- 1 ) 2 = —R 2 m 2 -{: d {d-l)+ (rf ~ 1)2 , (2.26) 

v ' ' fc=-(d-5)/2 
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where we have used eq. (2.9) in the second expression for b. 7 Further for d = 3 and 4, 
the product appearing as the last factor in f(y) should be omitted. 

Even though eq. (2.25) for ^-function is only valid for Re(z) < 1, this will be suffi- 
cient to compute the Renyi entropy in the present context since according to eqs. (2.4) 
and (2.17), we only need to know ^-function in the vicinity of z = 0. Note, however, 
that general expression which is valid for all values of z can be found in [31]. 

The ^-function approach is similar to dimensional regularization in that no power 
law divergences will appear with this method. Hence the leading contributions to 
the entanglement entropy of the form given in eq. (1.1) are somewhat obscure in this 
framework. 8 In contrast, evaluating the entanglement (and Renyi ) entropy using 
eq. (2.8), where 5 directly cuts off the UV end of the t integral, produces explicit power 
law divergences as appear in eq. (1.1). We illustrate these differences by applying both 
approaches in the examples below. Since the form of the heat kernels and ^-functions 
is different in even and odd dimensions, we consider these cases separately. 

Odd dimensions 

We start to implement eq. (2.8) in the special cases d = 1 and 3 to illustrate how the 
divergent 'area law' and subleading terms emerge, as well as the universal 'area' terms 
(1.3). We then consider the ^-function method (2.17) to evaluate the finite contribu- 
tions to the Renyi entropy for general odd d. 

d=l (D = 3): 

In this case one has to substitute eqs. (2.18) and (2.23) into eqs. (2.4) and (2.8) 



Note that no £ dependence appears in the expressions above because the curvature 
scalar vanishes on S 1 , i.e., rriQ = m for d = 1 in eq. (2.9). The divergent term is, of 
course, the expected 'area law' contribution. It originates from the n = summand in 
eq, (2.23) and thus is independent of the cross-section of the waveguide geometry, e.g., 
the same term arises in eq. (B.ll) for a hyperbolic waveguide. 

7 This notation is not ideal in the following where we focus on the limit mR 3> 1. Hence let us 
resolve the possible ambiguity by adding that b ~ i mR in this limit. 

8 Of course, this also illustrates the sensitivity of these contributions to the details of the UV 
regulator. 





(2.27) 
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In examining the finite contribution above, we first note that it does not have 
the simple form expected in eq. (1.3). However, we note that the full calculations in 
[1] resulted in a similarly complicated expression and the simple universal term only 
emerged in the large mass limit. Hence we consider the finite term above in the limit 
mR > 1, 

S$ inite = -±±£ (2,R m - 2e~ 2 ™ R - e~ 4 ™* + ••■). (2.28) 

Here we see that the leading term has precisely the form expected in eq. (1.3) with 
= 2ttR and D = 3. We emphasize that above expression describes the Renyi 
entropy and can be evaluated for any a. The entanglement entropy is recovered by 
substituting a = 1, in which case the pre-factor becomes —1/12 precisely matching the 
coefficient given in eq. (1.4) for D = 3. Here we also find higher order terms suppressed 
by exponentials exp(— 2nnmR). Similar exponential terms were found in [1] but the 
precise numerical prefactors do not agree for the waveguide geometry studied there and 
for the present cylindrical waveguide. Of course, these terms only become important 
when the Compton wavelength of the scalar is comparable to the size of the cross- 
section of the waveguide. Hence it seems these contributions are probing the topology 
of the background geometry. In fact, in eq. (B.12), we find that no such exponentials 
arise when the cross-section is H 1 ~ R 1 . 



d = 3 (D = 5): 

Taking the limit 9 — > in eq. (2.21), yields 



^3 } (t,0) 



D t/R 2 



{Ant) 3 / 2 



n 2 R 2 n 2 



1 - 2 



t 



(2.29) 



Combining this expression with eq. (2.18) to form K$ and in turn, substituting the 
result into eqs. (2.4) and (2.8) yields 



c(o) 



Vol(§ 3 ) ^ ^e-< m2+ ^ R2 )K^\t,0) 



(2.30) 



1 + a Vol(§ 3 



12 a (4tt) 3 / 2 



3<5 3 



2 6£ - 1 \ 1 
— I m + = I 



R 2 J 5 



a 



96 7r 2 a 



g {0) (2ir^J(mR) 2 + 6f - 1 



with Vol(§ 3 ) = 2tt 2 R 6 and 



g®{x) 



6' 



-x 



x logfl 



2Li 3 ( e - x ) -2xLi 2 (e" 



(2.3i; 



where Li 3 (z) and Li 2 (^) are the standard polylogarithms. In eq. (2.30), we see the 
expected area law contribution proportional to Vol(S 3 )/5 3 , as well as a subleading 
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divergences proportional to 1/8. Focusing on the finite contribution in the limit mR 3> 
1, we find 

cm _ i+« , ( m z 1 m - 1) m m - 1) 2 i i 

^finite ~ ^ Q ^ ^ + 2 r2 + g ^ ^ + • • ■ 

-- — e~ 2 ™ R +--- ) (2.32) 

where At, = Vol(S 3 ) = 2ir 2 R? is the area of the entangling surface. The leading term 
in this expansion has the form expected from eq. (1.3) and setting a = 1 to recover the 
entanglement entropy, the prefactor becomes l/(72n) which precisely matches the co- 
efficient given in eq. (1.4) with D = 5. There are two classes of sub leading terms: First 
there is an expansion in powers of l/(mR) 2 , which produces terms where the prefactor 
becomes ^4sm 3 / (mR) 2n . These contributions have precisely the form expected for the 
curvature terms described in eq. (1.5). Second, there are contributions with exponen- 
tial factors exp(— 27rn mR), similar to those found with d — 1. As explained above, it 
appears that these contributions probe the topology of the waveguide geometry. Note 
that, as illustrated in eq. (2.32), some of these exponential contributions contain odd 9 
powers of 1/R, which emphasizes that these terms cannot be given a simple geometric 
interpretation, as in eq. (1.5). 



General odd d > 3: 



As already noted in footnote 6, the logarithmic divergence in eq. (2.17) vanishes for 
odd d since C§° (0) = 0. Hence, evaluating the finite term in the Renyi entropy using 
eqs. (2.4), (2.17), (2.18) and (2.25) yields: 



o(0) 
a, finite 



-D^i±^Vol(§ rf ) 



d( 



(o) 

K d 



24a 



Vol(H d ) dz 



z=0 



Substituting eq. (B.4) for d5 above, then produces 



.2 d ~ 2 T{d/2) [°° f(ib-y)dy \ 

1 _ e 2w(y-ib) J 
(2.33) 



7T—R d 



a(0) 
a, finite 



1 1 + q Vol(S d ) / 1 ( V^ /2 jo) 2fc+1 sec(fc7r) 
j 24a w *? R d\2*-ir(d/2) h 2k + 1 

f(ib-y)dy 



k=0 



+i2 d - 2 T(d/2) 



I _ e 2ir(y-ib) 



3 Here, we mean odd powers of 1/R multiplying As, which itself contains a factor of R 3 



13 



where g^l = 0, #13 = 1 an d for odd d > 5 are defined by 

, o 2 (<*-5)/2 (rf-l)/2 

- 2 +(V) ] n (- 2 +^)= e ( 2 - 35 ) 

i=o fc=o 
From this expression, it is useful to note that 

A, = l and ^3 =^(d-l)(d-2)(d-3), (2.36) 

2 >" 2 >" 24 

as well as = 0. We may then use these expressions to expand S^ inite in the limit 



mR » 1 and in doing so, we find 

(0) 1 + q (-I)Vtt / D _ 2 (D-2) 2 (D-3)(6e-l)m^- 4 

^« 24a (47r) V r ( D/2 )^ E r + 12 ^ 

(2.37) 

where as before At, = Vol(S ) is the area of the entangling surface. Of course, setting 
D = 5, eq. (2.37) simply produces the first two terms in eq. (2.32). For general D, 
we may set a = 1 to recover the entanglement entropy and we see that the leading 
term above is the precisely the area term expected in eqs. (1.3) and (1.4). The next to 
leading term in eq. (2.37) introduces a new universal contribution to the entanglement 
entropy which matches the form shown in eq. (1.5) with n — 1. This contribution can 
be interpreted as 

su* = U " I) lD ' scalar ^ ^ nW mD ~ A (2 ' 38) 

where 71(h) is the Ricci scalar of the metric induced on the entangling surface and the 
coefficient jD^caiar is precisely that given in eq. (1.4). Note that we should only consider 
this term for odd D > 5. Of course, there are several other independent curvature terms 
which could in general contribute at this order but with an appropriate choice of basis, 
the extra terms all vanish for the waveguide geometry studied here — see section 4 for 
further discussion. 

Even dimensions 

Following our discussion of odd d, we first consider the special value d = 2 here and 
evaluate all UV divergences using eqs. (2.4) and (2.8). Then for general even d > 4, 
we apply the approach of ^-function regularization. This method eliminates power law 
divergences, while keeping the universal terms, i.e., logarithmic divergences, as well as 
finite contributions to the Renyi entropy. 
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d = 2 (D = 4): 



In the limit of coincident points, we get from eq. (2.24) 

1 oo 

K&\t) = — 2 £(2Z + l)e-V« . (2.39) 
i=o 

Applying the Euler-Maclaurin formula, i.e., 

E F «) - jT + F( ° )+ / (00) + £ ^ (F^)(oo) - F(^)(0)) 



(2.40) 



leads to the following expansion 



Combining this expression with eq. (2.18) yields as in eq. (2.11). Then substitut- 

ing the result into eqs. (2.4) and (2.8) yields 



S< 0) = ^Vol(§ 2 ) r^e-<^ + W*)K®{t) 
24a t 



1 + a 



(jn 2 + 6 ^ 3R2 1 ) \°g( mS ) + 



25 2 \ 3R 2 



(2.42) 



where As = Vol(§ 2 ) = 4nR 2 and ellipsis denotes finite terms. Upon setting a = 1, we 
recover the entanglement entropy and the first term is recognized as the standard 'area 
law' contribution. The logarithmic contribution proportional to m 2 matches the area 
term given in eq. (1.3) with D = 4. With £ = 1/6 and m 2 = 0, the theory (2.7) under 
consideration becomes a conformal scalar. In this case, one can verify that the loga- 
rithmic contribution above (which vanishes) matches the expected result from eq. (1.2) 
for a conformal scalar — see section 4 for further discussion. 

General even d > 4: 

According to eqs. (2.4) and (2.17) for even d, the universal term is proportional to 
Cg^(0). We can evaluate the latter using eq. (2.25). Further, we see that since sin(7rz) 
vanishes at the origin, the latter expression simplifies to 



® (o) = <- i)j/2 !S c » i(o) • (2 ' 43) 
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Hence, using eq. (B.6) for the scalar (^-function on M. d , we find that universal contribu- 
tion is given by 



S 



(o) 

a,univ 



1 + a 



-l)^Vol(S rf ) {d ^ 2 
12a (An) d / 2 T(d/2)R d f-^ 



h 



(0) 

k,d 



U2\k+l 



4 



oo j.2k+l 



d 2ttx 



-dx 



where h 



(o) 

0,4 



1/4, hi 



(o) 



- 1 and ftj^ for even d > 6 are defined by 



(d-5)/2 



(<f-2)/2 



n (- 2 +/)= e c 

J=l/2 



(0) x 2fc 



fc=0 



Given this definition, it is useful to note that 



h 



(0) 



1 and h 



(o) 



_(d-l)(d-2)(d-3). 

:(o) 



log(m5) , 
(2.44) 

(2.45) 
(2.46) 



Using the latter two expressions in an expansion of S^j inite in the limit mR » 1 yields 



S. 



(o) 



(-l)^Us 



12a (47r)Trr(iy2)( 



D _ 2 , (D-2) 2 (D-3)(6e-l)m 
m H 



D-4 



12 



i? 2 



f. . . J log(m5) . 



(2.47) 

As usual = Vol(S d ) is the area of the entangling surface. Note that setting D = 4 
(i.e., d=2), eq. (2.47) reproduces the universal term calculated above in eq. (2.42). 
Again with a — 1, the above reduces to the entanglement entropy and we see that the 
leading term is the precisely the area term expected in eqs. (1.3) and (1.4). Similar 
to the discussion for odd d, The next to leading term in eq. (2.47) introduces a new 
universal contribution to the entanglement entropy which again matches the form shown 
in eq. (1.5) with n = 1. We can write this contribution as 



S„ 



D-2 



S g ) ID, sea 



d D ~ 2 a Vh K(h) 



m 



D-4 



\og(m8) 



(2.48) 



where lZ(h) is the Ricci scalar of the metric induced on the entangling surface and the 
coefficient ^D,scaiar is given in eq. (1.4). Here we should only consider this term for even 
D > 4. Of course, this expression is reminiscent of eq. (2.38) for the case of odd d. 
If we evaluate the entire expression (2.44) for d = 4 (D — 6), we obtain 



s£L,(d = 4) 



1 + a Vol(S 4 ) 
a 192 7T 2 



m 



-m 4 + 2(6£-l)-^+(72£ 



24£+— I — 
^ 15 J i? 4 



\og(m5) . 
(2.49) 



This example illustrates that the curvature contributions in eq. (1.5) extend up to 
n = d/2 for even d, as can also be seen by directly examining eq. (2.44). In contrast 
to the case of odd d, these two equations also show that there are no exponentially 
suppressed terms in the universal contribution for even d. 
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2.2 Renyi entropy for a massive fermion 

In this section, we construct the partition function for a spin-| field living on the 
Euclidean manifold Ai a = C a x S d and use this result to evaluate the corresponding 
Renyi entropy. Our spinor notation is reviewed in Appendix A, whereas the action 
under consideration is given by 



S(i/j,i/j)= (i/jfip + mifjvp) . (2.50) 

In this case, the massless wave operator appearing in the heat kernel (2.8) is 

U-1/ ] = -1/\ (2.51) 

which we refer to as the 'iterated' Dirac operator, following [32]. Since the Dirac 
operator is a nondiagonal matrix, the separation of variables in eq. (2.11) is not obvious 
and so we prove that this equation still holds here in Appendix A. The argument there 
rests on the structure of the heat kernel for the iterated Dirac operator on § d , which 
we review next. 

Let us first consider the case of odd d and further for simplicity, let us assume 
that one of the points coincides with the north pole of E> d . In this case, the heat kernel 
reduces to [32] 

K^V^my)^-^^)'^-)- 1 £ (-1)"^. (2.52) 

N ' n=—oo v ' 

where y is an arbitrary point on the sphere. The angle of latitude for this point is 
designated as 9 and then 9 n = 9 + 2irn. Finally U(y) is the spinor matrix which 
parallel propagates a spinor from the given point y to the north pole. Similarly, in the 
case of even d, the heat kernel becomes [32] 

K^(t, y) = U(y) cos \ (^j^)' (cos \ )-'/n*. <^3) 

where 

* {<> ' t] = (i^W^m) „5L /. VcosO-cos^ 6 " # ■ (2 ' M) 
with n = + 27m. 

The structure of the spinor matrix U(y) can be found in [32] 10 but these details will 
be not important here because we are only interested in the limit of coincident points. 



10 



We also refer the interested reader to eq. (A. 12). 
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In this limit, U (y) simply reduces to an identity matrix. We might note that if we were 
considering the heat kernel on S d alone, the dimension of this identity matrix would be 
2 L<^/ 2 J ^ ^ e ; the dimension of Dirac spinors in d dimensions. Of course, here S d is part 
of the larger manifold Ai a and so the dimension of U(y) is actually 2L D / 2 J. However, 
following the conventions introduced in Appendix A, we treat the spinor trace on the 
right-hand side of eq. (2.12) as though we separately tracing over the spinor spaces of 
the two component manifolds. That is, we calculate the two spinor heat kernels on C a 
and S d separately and then simply take their product in eq. (2.12). 

The spinor heat kernel on the cone C a is readily evaluated as [29, 30] 

Tr< /2) (t) = ~ (l - « 2 ) +aTrKg /2 \t) . (2.55) 

As noted above, this result accounts for the trace over the two-dimensional spinor 
indices on the cone C a . 

Unfortunately, in the present case, we are unable to apply the ^-function approach, 
which would have allowed a systematic evaluation of the universal contributions to the 
Renyi entropy for general d. In particular, while the spinor ^-function is known for 
the hyperbolic space H d [32], the spin-- counterpart of eq. (2.25) relating these to the 
desired (^-functions on § d is unavailable. Hence, in the following, we limit ourselves to 
considering a few special cases, i.e., d = 1, 2 and 3. In each case, the Renyi entropy 
is determined by simply substituting eq. (2.52) or (2.53), along with eq. (2.55) into 
eqs. (2.4) and (2.8). Given these expressions above, the generalization of the following 
results to higher dimensions would be straightforward. 

Let us add that we consider fermions on the hyperbolic waveguide C a x W d in Ap- 
pendix B. In this case, we can follow the (- function approach using the results of [32]. 
This allowed us to produce results for the Renyi entropy of spin-| fields for general d 
and in particular, gave the general coefficients appearing in eqs. (1.6) and (1.9). 

d=l (D = 3): 

In this case using eq. (2.52), we find 

n=— oo 

^g-W-l-)]). (2.56) 
In the large mass limit (i.e., mR 3> 1), the finite contribution becomes 



S2£L = ~- (2*R m + 2e" 2 ™* - e~^ R + ■■■). (2.57) 
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Hence the leading term has precisely the form expected in eq. (1.3) with At, = 2ttR 
and D = 3. The entanglement entropy is recovered by substituting a = 1, in which 
case the pre-factor matches that given in eq. (1.6) for D = 3. Note that this area term 
is accompanied by higher order exponential terms, similar to those found in eq. (2.28). 

d = 2 (D = A): 

In the limit of coincident points, eq. (2.53) yields 

D +00 „2L 2 2 

< /2) (0 = A 1/2) (0,t) = £ jj^ + nn) cot0 e-^^# . (2.58) 

We may note that the integrand is regular at the lower bound. While cot <fi has a 

simple pole at (ft = 0, the coefficient of this pole is proportional to Xln^oo ne_ * 
which vanishes because the summand is odd in n. Redefining the integration variable 
</>—>•</>/ \/t, we have 

D +oo c-Z-~ / \2 

K &'» = 2-5^175 E / (#+ 9 cot(v^)e- R "^ # . (2.59) 

n=— oo ^ * 

From this expression, it is obvious that the UV divergences in eq. (2.8) come entirely 
from the n = term. Any terms with n ^ contain an exponential factor e ~' K ' lR2n ' 1 l t 
which smoothes out any potential singularities at t = 5 2 . Therefore, to extract the 
structure of the UV divergences, and in particular, the log(mJ) contribution, it is 
enough to examine only the n = term, i.e., 

*S /2) (*)= o stLiTS /^0cot(^)e-^ 2 #. (2.60) 



2 7T 3 / 2 tV2 , 

The remaining terms only contribute to the finite part of the Renyi entropy. 

Now we expand the integrand in eq. (2.60) in the vicinity of t = and find that 
it is sufficient to keep only first two terms since the rest do not lead to singularities in 
the limit 5 -»■ 

-(1/2) /.n 1 f°° ( , 1 ji * -~ 2 . 



i^(t) = — / (l -x 2 + ...)e- x dx + ... = - + ...(2.61) 

s v ; 2tt 3 /H J K 3R 2 J Ant 24nR 2 v ; 

Combining this result with eq. (2.55) to form iirj^ and substituting the result into 
eqs. (2.4) and (2.8), we find 

Sr» = ^Vol(S 2 )^fe-'™ 2 Jf<;' 2 »( ( ) 

= ^^(i + ( 2m2+ 3^) l0 ^ + -) < 2 ' 62 ' 
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where Ay = Vol(§ 2 ) = 4ttR 2 and ellipsis denotes finite terms. For a = 1 the first term 
represents the standard 'area law' in the entanglement entropy, whereas the next term 
proportional to m 2 is precisely the D = 4 case of eqs. (1.3) and (1.6). Similarly, the 
term proportional to 1/R 2 matches eqs. (1.7) and (1.9) with D = 4. 

d = 3 (D = 5): 



Taking the limit 9 — > with d = 3 in eq. (2.52), yields 



1 00 2 / 



2R 2 



(2.63) 



Combining the above with eq. (2.55) to form K^J^ and substituting the result into 
eqs. (2.4) and (2.8), yields 

1 + a 



Q(l/2) 



Vol(§ d 



dt 



12a v ' Jp t 
1 + a Vol(§ 3 ) r 1 



6 a (4vr) 3 / 2 



36< 



KM 



- 77T + 



2i? 2 /5 



where 



1 + a 
48 n 2 a 



g {1/2) (27rRm) , (2.64) 



g (i/2) {x) = * + -« + n 2 log{1 + e - B) + Li3 (_ e -x) + xL i 2 (-e" 
6 2 

Expanding the finite contribution above in the limit mR 3> 1, we find 



7T 



(2.65) 



? (0) 
a, finite 



a 



a 



Ay m 



3 m 



3 m 



-2wmR 



4 i? 2 2vr 2 R 2 



(2.66) 



where = Vol(§ 3 ). Setting a = 1 to recover the entanglement entropy, the pre-factor 
becomes l/(367r). In this case, the leading term in this expansion is the area term 
(1.3) with precisely the coefficient given in eq. (1.6) for D = 5. Also, the next term 
matches eqs. (1.7) with the coefficient given by eq. (1.9) for D = 5. Further, the above 
expansion also reveals contributions with exponential factors exp(— 2irn mR), similar 
to those found previously for odd dimensions. 



3 A calculation at strong coupling 

In this section, we are going to use gauge/gravity duality to study the universal 'area' 
contribution to the entanglement entropy of M = 2* gauge theory [20, 21] at strong 
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coupling. The latter is a massive deformation of the four- dimensional M = 4 super- 
Yang-Mills (SYM) theory, which is commonly studied in the AdS/CFT correspondence. 
So let us begin by giving the field theoretic description of the mass terms which appear 
in this context. The M = 4 SYM theory includes a gauge field A^, four Majorana 
fermions if) a and three complex scalars fa, all of which are in the adjoint representation 
of the U(N) gauge group. Now there are two independent 'mass' terms which can be 
used to deform the SYM theory [24] 

5C = -2 J d 4 x [m 2 b 2 + m f 3 ] (3.1) 

where 

O 2 = iTr(|0 1 | 2 + |0 2 | 2 -2|0 3 | 2 ) , (3.2) 

O s = -Tr (i ^ 2 - V2g YM 3 [0i, 4>\] + V2g YM <f) 3 [4, 2 ] + h.c.) (3.3) 

+ jjm / Tr(|0 1 | 2 + |0 2 | 2 + |0 3 | 2 ) • 

The M = 2* gauge theory results when we set mj = rrif and the effect is to give 
mass to the M = 2 hypermultiplet comprised of 0i i2 and ^1,2 • At this point, we note 
that the dimension-two operator (9 2 contains an unstable mass term for the scalar 3 
- a typical characteristic of such superconformal primary operators. This negative 
mass-squared contribution from O2 is precisely canceled by the positive contribution 
in 3 when mj = rrif and hence <p 3 is left massless in the supersymmetric theory. 
In the following, we will not restrict our attention to the supersymmetric theory and 
instead we set rrib and rrif to independent values. In this case, one may then worry 
that the resulting theory is unstable, i.e., when > rrif. Determining the end-point 
of this instability is the question of understanding the infrared behaviour of the RG 
flow induced by the mass terms. 11 However, recall that the universal contribution to 
the entanglement entropy appears with a logarithmic dependence on the UV cut-off, 
log S, as shown in eq. (1.3), since we are studying a four-dimensional gauge theory here. 
Furthermore, as demonstrated in [2] and as we will explicitly see below, the coefficient 
of this logarithmic term is only determined by the UV properties of the RG flow and is 
completely insensitive to the IR details. Hence, any such instability is of no consequence 
to the following calculations. 

n The IR theory can be stabilized by introducing a finite temperature, e.g., [22, 23]. However, such 
a modification of the IR state would again not modify the coefficient of the logarithmic contribution 
in the entanglement entropy [2]. 
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The dual holographic theory consists of five-dimensional Einstein gravity coupled 
to a pair of scalars, a and x, as described by the following action [20, 21]: 

h = [ <U*yT9 £5 = tttV / d ^^9 i R ~ 3 ( 9a f - (d X ) 2 - V(a, X )] , (3.4) 

J Mb l07TGr 5 J M5 

where the potential takes the form 

V(a, X ) = -^e" 2a - cosh x + ± e *> sinh 2 x . (3.5) 

Given the above expression, it is trivial to show V(a = 0, x — 0) = ~~ 12/ L 2 and so we 
see that L corresponds to the curvature scale of the AdSs vacuum solution. As it will 
be needed below, we also note that in the present conventions, Newton's constant is 
given by 

< 3 - 6 > 

where N is the rank of the U (N) gauge group in the boundary theory. As is well-known, 
the AdS/CFT correspondence relates the asymptotic boundary behaviour of the scalars 
to the couplings and expectation values of the dual operators in the boundary theory. 
Here, a is dual to the "bosonid mass coupling ml and the corresponding operator (D 2 , 
given in eq. (3.2). Similarly, x is dual to the 'fermionid mass rrif and the operator O3 
in eq. (3.3). 

In this holographic context, we follow the now standard approach to calculating 
entanglement entropy [8, 16]. That is, given a spatial region V in the boundary theory, 
the entanglement entropy between this region and its complement is given by the 
following expression evaluated in the bulk spacetime: 

S(V) = ext [A(v)} . (3.7) 

Here v is a (three-dimensional) bulk surface extending out to asymptotic infinity such 
that its asymptotic boundary dv matches the 'entangling surface' dV in the boundary 
geometry. The symbol 'ext' indicates that one should extremize the area over all such 
surfaces v. 12 

Following [24] (see also [22, 23]), we adopt the following ansatz for the bulk solution 

ds 2 = % \e 2A ^(-B(z) 2 dt 2 + dx 2 ) + dz 2 \ (3.8) 
a = at{z), X = x(z). 



12 If eq. (3.7) is calculated in a Minkowski signature background, the extremal area is only a saddle 
point. However, if one first Wick rotates to Euclidean signature, the extremal surface will yield the 
minimal area. 
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This solution reduces to simply AdSs with A — — B — a — x- 111 general, turning on 
the bulk scalars induces a holographic RG flow which is encoded in the metric function 
A(z). The function B(z) provides a potential blackening factor in the time compo- 
nent of the metric, which also allows the above ansatz to describe finite temperature 
situations. 

In order to identify the log 5 term in the entanglement entropy, we will only need 
to consider the asymptotic or small z region of the bulk solution, which describes the 
UV behaviour of the dual gauge theory. The desired asymptotic solution was found in 
[24] and takes the form 

a(z) = z 2 (a ,i log(» + a ,o) + 0(z 4 log 2 z) , 

X(z) = z xo,o + z 3 Qxo.o log z + X 2,o^j + 0(z 5 log 2 z) , (3.9) 

A(z) = -^z 2 + 0(z 4 log 2 z), 
B(z) = 1 + B ifi z 4 + 0(z 6 ). 

Of particular importance, the coefficients, cto,i an d Xo,0; correspond to the field theory 
masses [24] 13 

2 

ao,i = -r and Xo,o = ^f- (3.10) 
o 2 

Similarly, the coefficients ao,o and X2,o ar e related to the expectation values of O2 and 

O3, respectively, while -B 4)0 yields the energy density. We might add that the exact 

solution describing the supersymmetric flow (with rrif — m& = m) is known [20] and 

the corresponding asymptotic expansion then becomes [22, 23], 

«(*) = (log(^) + I) + 0(mV) , 

X (z) = ^* + ^ 3 (log(^) + i) + 0{m 4 z 4 ) , (3.11) 

A{z) = -^m 2 z 2 + 0(m 4 z 4 ), B(z)=0. 

Turning now to the entanglement entropy, our goal is to compare our holographic 
results to those found at weak coupling. To facilitate this comparison, we must perform 
the holographic calculation for a 'waveguide' geometry analogous to those considered 
in the previous section. Given the above ansatz (3.8), the boundary geometry is just 
flat space, but if two of the spatial coordinates are periodic, this geometry can be 

13 With respect to conventions in [24], we have already extracted the L dependence in the metric. 
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interpreted as a toroidal waveguide M 2 x T 2 . That is, we will identify the spatial 
coordinates, x l and x 2 , with a period H where H ^> l/m^, 1/irif. The entangling 
surface £ is then chosen as x 3 = (and t = 0). According to eq. (3.7), we must find 
the extremal bulk surface that connects to this entangling surface as z — > 0. However, 
the high degree of symmetry here dictates that the extremal surface will simply fall 
straight into the bulk geometry. That is, the desired bulk surface is simply given by 
v = {x 3 = , t = 0}. Next we must evaluate the area of this surface, however, as usual, 
the latter must be regulated by cutting off the integral at the UV regulator surface 
z = 5. We are particularly interested in identifying any contributions proportional to 
log S and so it is sufficient to consider the asymptotic solution (3.9) in calculating the 
area. 

A(v) = J dx 1 dx 2 dz y/g^g 22 9zz = H 2 Jdz^ e 2A{z) 

"^/^ (i-^o^ + O^log 2 ^ (3.12) 
- ^ + \h 2 L 3 x 2 0fi log 5 + 0(5 2 log 2 5) . 

Now we may combine this result with eqs. (3.6) and (3.10) to write the entanglement 
entropy as 

A(v) N 2 Av N 2 A S 2 . , „ , 1 , 

where At, = H 2 is the area of the entangling surface in the boundary theory. Of course, 
the leading contribution here is the expected 'area law' contribution, which does not 
yield any universal information. The next term yields the desired logarithmic term, 
with the same general form (1.3) as found from free field calculations. We may observe 
that, as is already evident in eq. (3.12), neither of these UV divergent terms depends 
on the higher order coefficients in the asymptotic expansion (3.9). Hence as expected 
[2], our calculation explicitly shows that these contributions to S* EE are insensitive to 
the IR details of the holographic RG flow. Of course, we may also observe that only 
the fermionic mass mj appears in the logarithmic term and the result does not depend 
on the bosonic mass mj,. 

While the above result applies to the gauge theory at strong coupling, we can 
calculate the analogous contribution to the entanglement entropy in the weak coupling 
limit. In this case, we simply consider the free field limit and apply the results given 
in eqs. (1.4) and (1.6) in eq. (1.3). If we first turn to (9 2 in eq. (3.2), we see that there 
are three complex scalars which acquire masses. Since each of these fields is in the 
adjoint representation of the U(N) gauge group, each fa effectively contributes 2iV 2 
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real scalars in the free field limit. Hence using eq. (1.4), the total coefficient in the 
logarithmic contribution becomes 

1 99 N 2 fml ml 2ml\ 

74 b = > 2N 2 m 2 = — + — =0. 3.14 

' ' 24tt ^ 1 12tt V 3 3 3 J K J 

Therefore in the weak coupling limit, the bosonic mass term does not in fact contribute 
to this universal area term in the entanglement entropy. Remarkably, this agrees with 
our strong coupling result (3.13), in that the bosonic mass m& does not appear in this 
expression. 

Next, we consider the contribution of O3 in the free field limit. In this case, all three 
scalars again acquire masses and the two fermions tpi and ip2 also acquire a mass. In the 
massless limit, ip± and ip2 are two independent Weyl fermions but here they combine 
as a single massive Dirac fermion. Again the latter are in the adjoint representation 
and so effectively we have N 2 Dirac fermions in the free field limit. Hence combining 
the results in eqs. (1.4) and (1.6), the total coefficient in the logarithmic contribution 
becomes 

1 2 2 1 2 f 2m2 f 2?Ti f 2 ™f\ N 2 2 , 

Comparing the above with the strong coupling result (3.13), we see that this coefficient 
is larger by a factor of three than that appearing at strong coupling. We can combine 
eqs. (3.14) and (3.15) to evaluate the coefficient in the supersymmetric theory with 
rrif = mj = m: 

N 2 

l±,susy = J4,f{ m f = m)+ 7 4 ,s(m 6 = m) = — m 2 . (3.16) 

47T 

Similarly, this choice of masses does not effect the holographic result (3.13) and so the 
discrepancy observed above between the coefficients in the strong and weak coupling 
limits extends to the supersymmetric case. 

Note that the strong coupling coefficient seems to match the contribution of the 
fermions alone at weak coupling. That is, we would have found agreement between 
the two limits if the 'fermionic' mass term C 3 only gave mass to the fermions ipi and 
'ip2- However, as is evident from eq. (3.3), this dimension-three operator also contains 
a mass term for all three scalars. The latter is somewhat unusual as it is a 'coupling- 
dependent' correction induced at finite mass, i.e., mjOs contains a contribution of 
order m 2 . However, the presence of the additional interactions in eq. (3.3) is dictated by 
the supersymmetry algebra and the global SO(6)r symmetry of the M = 4 theory [24]. 
Further, one can directly detect these scalar masses [34] by uplifting the asymptotic 
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solution (3.9) to ten dimensions [22] and then examining the potential felt by a probe 
D3-brane, e.g., [21]. Of course, we should also note that beyond the fermion and scalar 
mass terms, 0% also contains trilinear couplings between the hypermultiplet scalars <pi^ 
and the gauge multiplet scalar </> 3 . Hence m,f does not simply parameterize the masses 
of various fields but also plays a role as the coupling of a new cubic potential term in 
the interacting theory. The latter may well be the source of the discrepancy observed 
above. 

4 Discussion 

In this paper, we studied the entanglement entropy entropy for a variety of field theories 
on waveguide geometries with a spherical or a hyperbolic cross-section. Our calculations 
confirmed the appearance of a universal area term (1.3), which was first identified in 
[1], as well as reproducing the precise coefficient (1.4) for a free massive scalar field. 
Our analysis of the scalar fields in section 2.1 also included a nonminimal curvature 
coupling ^C,R(p 2 and we found the previous area term is not affected by this new 
coupling. In section 2.2, we also considered a free massive fermion and we identified 
the same universal area terms as in eq. (1.3) with the coefficient given in eq. (1.6). This 
reproduces a result given previously in [19]. 

Curvature contributions: 

By considering waveguide geometries with a curved cross-section, we were also able 
to consider new curvature contributions in the entanglement entropy of the schematic 
form shown in eq. (1.5). Such terms were first identified with holographic techniques [2], 
however, we can now present a clearer understanding of the origin of such contributions 
to entanglement entropy, following the perspective given in, e.g., [35-37]. As noted in 
the introduction, the calculation of entanglement (or Renyi) entropy in QFT generically 
yields a UV divergent answer because the result is dominated by the short distance 
correlations in the vicinity of the entangling surface E. Hence the calculation must be 
regulated by introducing a short distance cut-off 5, and the result typically contains 
a series of power law divergences, as shown in eq. (1.1) for a QFT in D spacetime 
dimensions. 14 Of course, the leading contribtuion has a geometric structure, in that it 
corresponds to the famous ' law' term [4] with: 




(4.1) 



Of course, with even D, the cd term may appear with a logarithmic divergence. 
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where ^ a b is the induced metric on the entangling surface. However, if we are working 
with a covariant regulator (in a relativistic QFT) and assuming the short-distance cut- 
off is much smaller than any scale defined by the couplings of the QFT, i.e., <C 1, 
then in fact all of the coefficients of divergent terms in this expansion (1.1) exhibit a 
similar geometric structure. For example, the second coefficient may be written as 



where, e.g., ^(7) denotes to the intrinsic Ricci scalar of the entangling surface. 15 This 
geometric character of the coefficients naturally follows from the fact that the UV 
divergences are all local. 

The dimensionless coefficients d 2 k,a above will of course depend on the detailed 
structure of the underlying QFT. However, because the d 2 k, a are dimensionless, we can 
write their dependence on any mass scale fii in the QFT in terms of the dimension- 
less combination faS, i.e., d 2 k, a = ^2fc,a(/^)- Unfortunately, the coefficients appearing 
in the expansion above are scheme dependent. Clearly, if we shift 5 —> a5, we find 
d2k, a d 2 k, a = ot 2k ~ D d2k,a{otiii8) . Hence the regulator dependence here comes both 
from the implicit dependence on mass scales in the QFT and the 'classical' engineer- 
ing dimension of the individual coefficients c 2 k- Of course, the latter reasoning can be 
evaded in certain special circumstances. One well-known example, which was already 
noted in the introduction is the case of a CFT in an even number of spacetime dimen- 
sions [6-10]. Given the underlying field theory is conformal, there are no intrinsic mass 
scales and in an even spacetime dimension, the cp term still accompanies a logarithmic 
divergence. Hence the corresponding coefficients d^^a evade the above scaling argument 
and provide universal information about the underlying field theory. In fact, as illus- 
trated in eq. (1.2), these coefficients are proportional to the various central charges of 



Another example, which we might consider, is where the underlying field theory 
describes some renormalization group flow beginning at a UV fixed point. The theory 
that is probed by the UV singularities (1.1) in the entanglement entropy is a CFT 
perturbed by some relevant operators. Now let us consider the special case where the 
RG flow in the UV is controlled by a single 16 mass scale [i. By assumption fi5 <C 1 and 
so we can express the coefficients in terms of a Taylor series: d 2 k, a = ^2 n =o^2ka (^) n - 
At this point, it is straightforward to extend the previous discussion to show when the 

15 We refer the interested reader to [11] for a full explanation of the notation. 

16 With more than one scale, one could carry out the following analysis with the largest scale /xo 
while allowing the coefficients d^ a to be functions of the ratios [iij 




(4.2) 



+d, A Kl b Kl a + d^Kl a Kl b ] 



the CFT. 
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coefficient (%. a with n = D — 2k does not scale with a. Hence, we can expect that this 
coefficient provides universal information about the underlying RG flow. Ultimately, 
we may repackage this discussion to see that we have identified possible universal con- 
tributions to the entanglement entropy which take the form 

See ^ 4 D k - 2k) fi D ~ 2k [ d D ~ 2 a^ [("curvature")'- 1 } , (4.3) 

JT, 

which, of course, matches the form of the expressions given previously in eq. (1.5). 

Of course, the previous discussion can not be complete as we know that univer- 
sal contributions can also appear with logarithmic factors \og(fi5), as illustrated in 
eqs. (1.3) and (1.7). So in general, we must allow for such \og(fi5) factors in the ex- 
pansion of d,2k, a to account for these situations. At this point, we would like to note 
that holographic techniques were used to study the effect on entanglement entropy of 
perturbing a strongly coupled CFT with relevant operators in [2]. One of the inter- 
esting results there was to show that logarithmic contributions can appear in either 
even or odd spacetime dimensions, in situations where the operators had large anoma- 
lous dimensions. For example, perturbing by an operator with conformal dimension 
A = (D + 2)/2 would generate such a log 5 contribution in general. Further the 
holographic calculations in [2] show that such large anomalous dimensions may also 
introduce unusual powers in the expansion of the coefficients d,2k,a- 

The appearance of universal contributions as in eq. (1.5) were first uncovered with 
a holographic approach in [2] and with sufficient effort, these holographic calculations 
allowed their precise geometric form to be identified in a straightforward way. How- 
ever, this identification is much more difficult in the context of field theory calculations 
presented in this paper. In fact, the geometric nature of various universal contributions 
is easily confirmed as follows. First, we observe that for sufficiently large masses, our 
results are naturally be presented in terms of an expansion in powers of l/(mR) 2 , as 
could be anticipated from previous calculations [19, 38, 39]. Of course, this is pre- 
cisely in agreement with the geometric expansion discussed above or given in eq. (1.5), 
where the powers of 1/R 2n correspond to factors of ["curvature"]™. This identifica- 
tion is further confirmed by repeating the calculations for hyperbolic waveguides in 
Appendix B. There we found that the same expansion is produced up to the replace- 
ment l/(mR) 2 — > —l/(mR) 2 , as compared to the spherical waveguides in section 2. 
Of course, this sign corresponds precisely to the change in the sign of the curvature 
between these two families of geometries. 

While in general extracting the precise geometric form of these contributions is 
difficult with the heat kernel approach used here, the simplicity of our background 
geometries allowed us to identify the first such universal contribution. As shown in 
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eq. (1.7), it is simply an integral over the entangling surface of the Ricci scalar evaluated 
on this surface. In general, one expects that five independent curvature terms would 
contribute at this order, as shown in eq. (4.2). However, we have chosen a particular 
'basis' for these terms there such that the expression identified in eq. (1.7) is the only 
nonvanishing term for the waveguide geometries studied here. To be precise, the last 
two possible contributions in eq. (4.2) involve an integral of terms quadratic in the 
extrinsic curvature of the entangling surface. However in the present construction, the 
extrinsic curvature is precisely zero and so these contributions vanish here. Similarly, 
the second and third contributions in eq. (4.2) involve the background curvature in the 
space transverse to the entangling surface. In the present case, this transverse geometry 
is simply M 2 and so these contributions are also zero. 

The coefficient of the curvature contribution identified in eq. (1.7) is given in 
eqs. (1.8) and (1.9) for the free scalar and fermion fields, respectively. In the scalar case 
(1.8), we see that this coefficient depends on the non-minimal coupling £. Focussing 
on this linear £ dependence, we note that the coefficient is precisely that which would 
appear if we replaced the mass m in the corresponding area term (1.3) by the effective 
mass m appearing in the action (2.7), i.e., 

m °- 2 = (m 2 + £71 ) (D ~ 2)/2 = m D ~ 2 + ^^m D " 4 £ K + ■ • • . (4.4) 

Here we should keep in mind that these results, i.e., eqs. (2.38) and (2.48), were derived 
from expressions where we had assumed mR ^> 1. Hence it is reasonable to consider 
the Taylor series expansion above. It seems that this simple dependence originates with 
the appearance of m 2 , in the definition (2.26) of the parameter b 2 which characterizes 
the mass dependence of our heat kernels. We might also note that D = 4 seems 
distinguished in eq. (1.8) or (2.48) in that this contribution vanishes when £ takes the 
value for a conformal scalar. 

As an aside, we point out here that there are finite contributions to the entangle- 
ment (and Renyi) entropy which do not fall into the class of terms discussed above. 
For example, as shown in eqs. (2.28) or (2.32), our results for the spherical waveguides 
revealed certain 'topological' contributions which are exponentially suppressed. We re- 
ferred to these terms as topological because there were no analogous contributions for 
waveguides with a hyperbolic cross-section. However, we also see in, e.g., eq. (2.32), 
that there are contributions with inverse powers of the mass. Recently, it has been 
shown that such terms with negative powers of m also have a universal character, as 
they can be related to the universal contributions appearing in entanglement entropy 
in higher dimensions [19, 38]. Note that if we examine the universal log 5 contribution 
in even D, none of the individual terms in this coefficient have the topological character 
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or inverse powers of the mass described here. The same is true for the coefficients C2k 
of the power law divergences in eq. (1.1). 

Renyi entropy: 

Our calculations also extend the initial work of [1] by presenting expressions for the 
Renyi entropy, as well as the entanglement entropy. In general, our results for the 
Renyi entropy have the same structure as the entanglement entropy. In fact, for any 
of the geometries or field theories which we studied here, we can write 

1 + a 

S a = ^—S EE . (4.5) 
2a 

Hence in general, there are a variety of universal contributions in S a which take the 
same geometric form as described above or in eq. (1.5). Unfortunately, with the relation 
in eq. (4.5), the Renyi entropy would not provide any new information that is not 
already available in the entanglement entropy. However, it seems that this simple 
'factorization' of the Renyi entropy must be an artifact of the simplicity of both the 
background geometries and the QFT's studied here. Typically, the Renyi entropy has 
a more complex dependence on the index a than appears in eq. (4.5) as can be seen, 
e.g., in the holographic results of [40] or the results for disjoint intervals in [41]. Both 
of these examples also demonstrate that the Renyi entropy often contains far more 
information about the underlying field theory than the entanglement entropy alone. 

The form of the Renyi entropy in eq. (4.5) is reminiscent of well-known results 
derived in two dimensions [25]. We can compare our expressions for the Renyi entropy 
for D — 2 (d — 0) to those given in [25] as a check of our calculations. In that reference, 
the authors consider a two-dimensional CFT perturbed by a relevant operator which 
introduces a correlation length 1/m. Combining various expressions appearing there, 
the Renyi entropy becomes 

S = ~^ c log(m5) , (4.6) 
12 a 

where c is the central charge of the CFT defining the UV fixed point. In our approach, 
the result for a massive scalar can be derived by substituting eq. (2.18) directly into 
eq. (2.8) and the answer then takes precisely the form given above with c = 1 as is 
appropriate for a free scalar field (in the conventions of [25]). Similarly substituting 
eq. (2.55) into eq. (2.8) yields a result for a massive Dirac fermion, which again has 
the form given in eq. (4.6) with c = 1. Here, we might note that eq. (2.8) has an extra 
minus sign in the case of the fermion relative to the scalar. However, eq. (2.55) also 
has an extra minus sign in comparison to its scalar counterpart and so we obtained the 
desired match. 
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Strong coupling: 

Previously, ref. [2] used holographic techniques to study the effect on entanglement 
entropy of perturbing a strongly coupled CFT with relevant operators. In the gravity 
description, such an operator is dual to a scalar field and the holographic entanglement 
entropy (3.7) is modified by the backreaction of the scalar on the geometry through 
Einstein's equations. As noted above, one of the interesting results was that logarithmic 
contributions could appear in either even or odd spacetime dimensions, in situations 
where the operators had large anomalous dimensions. For example, perturbing by 
an operator with conformal dimension A = (D + 2)/2 would generate such a log 5 
contribution in general. Another result, which seemed to present a small puzzle, was 
that an operator with the dimension of a scalar mass term, i.e., A = D — 2, generated 
a log 6 term in even dimensions but only for D > 6. The puzzle was then that at weak 
coupling, i.e., with free field theories, such a logarithmic term appears for a massive 
scalar [1], as shown in eq. (1.3). 

Here in section 3, we considered a particular holographic framework where the 
AdS/CFT dictionary is very well understood, namely the four-dimensional M = 2* 
gauge theory. In this case, two mass operators with A = 2 and 3, as well as their 
holographic description in the dual gravity theory, are known exactly. The detailed 
description of the bosonic mass term (3.2) in the boundary theory seems to resolve 
the puzzle noted above. In particular, the intuition provided by weak coupling is that 
the coefficient of the log S term is proportional to m i m a case where several scalars 
acquire masses. However, this sum is precisely zero for the present mass term given 
in eq. (3.2). Hence rather than a puzzle, we have precise agreement between the weak 
and strong coupling results in this specific case. 

In examining the effect of the 'fermionic' mass operator (3.3) on the entanglement 
entropy, we found a contribution in eq. (3.13) of the general form A^mj log(m/<5), as 
expected an analogous weak coupling calculation. Unfortunately, if we compare the 
precise coefficient found at strong coupling with that appearing in the free field limit, 
there is a discrepancy by a factor of three. Note that the latter adds contributions from 
both the fermions and scalars which acquire masses when this operator is introduced. 
As shown in eq. (3.3), O3 also contains a cubic interaction between the scalar fields, 
which vanishes in the weak coupling limit, i.e., g YM — > 0. However, in the present 
context then, rrif does not only parameterize the masses of various fields but it also 
plays a role as the coupling of a new cubic potential term in the interacting theory. Of 
course, it is tempting to argue that the latter is the source of the discrepancy observed 
between the strong and weak coupling results. An interesting extension of the present 
results then would be to calculate the effect of the cubic interaction on the entanglement 
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entropy perturbatively when the gauge coupling is small but finite. 

We would like to contrast the above discrepancy with the recent results in [42]. 
There the author found in perturbative calculations, that the effect of the interac- 
tions on the entanglement entropy was to properly renormalize the mass appearing in 
eq. (1.3) so that it corresponded to the physical mass. Beyond the usual difficulties that 
one would encounter in extending such a calculation to strong coupling, our holographic 
calculation points out another difficulty in taking this limit. Namely, our strongly cou- 
pled boundary theory has no simple (quasi)particle excitations and hence we could not 
identify the 'renormalized mass' from a pole in a two-point function. Clearly, a better 
understanding is needed to appreciate the precise sense in which the entanglement en- 
tropy contributions identified in [1] are universal or alternatively, to unravel the precise 
information that these terms carry about the underlying field theory. It would be useful 
study further holographic examples where the precise definition of the relevant opera- 
tors is known in both the bulk gravity and boundary field theories. Another example, 
which would be interesting for this purpose, is the Cvetic-Gibbons-Lu-Pope solution 
[43], which describes an RG flow from a three-dimensional CFT [44]) in the UV to a 
gapped theory in the IR. 

Comparison with CFT's: 

As noted previously, certain universal contributions to the entanglement entropy are 
well known in the case of CFT's in an even number of spacetime dimensions [6-10]. 
While our calculations focussed on the universal contributions appearing with masses, 
they should also yield the expected CFT results in the appropriate limits. Hence it is 
interesting to compare our expressions with the expected CFT results as a check of our 
calculations. 

The universal contribution for a four-dimensional CFT [9] is given in eq. (1.2). 
Setting D = 4 in section 2, the corresponding entangling surface is E = S 2 . In our four- 
dimensional waveguide geometry, the extrinsic curvatures vanish, TZy, = 7£§2 = 2/R 2 



Moreover, the central charges for a massless conformal scalar (£ = 1/6) and a massless 
fermion are [45]: 




ftg 2 /3 = 2/(3R 2 ). Hence, eq. (1.2) yields 




(4.7) 




for s = , 
for s = 1/2 , 



and 




for s = , 
for s = 1/2. 



(4.8) 
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Therefore, we finally find 

S univ (s = 0) = , and S univ (s = 1/2) = -L log(5/R) . (4.9) 

lo 

Now we may compare these results with those derived in section 2. In particular, if 
we set m = and £ = 1/6 (and a = 1) in eq. (2.42), we see that the log 5 term in 
the entanglement entropy vanishes for a conformal scalar, in agreement with the above 
result. Similarly if we set m = (and a = 1) in eq. (2.62), we recover precisely the 
above expression for the universal contribution of a massless Dirac fermion. 

Our waveguide geometry also lends itself to using the approach of [10] to determine 
the universal contribution for a CFT in any number of dimensions. The only restriction 
of this latter approach is that the background geometry must have a rotational symme- 
try in the transverse space around the entangling surface, which is certainly satisfied 
in the present case. The result for the six-dimensional waveguide M 2 x S 4 is given in 
[11] as 

Si = ^§ i^A - 175! + 52B 2 - 592B-A log(mS) (4.10) 

where A, B\, B 2 , and B% are the four central charges which appear in the trace anomaly 
of a general CFT in D = 6. These central charges were found in [46] for a massless 
conformal scalar (£ = 1/5): 

5 28 5 2 

A = r , B] = — — ~ — - , Bo = ; — — r , and Br, = — - — — — ; (4.11) 

3-7!' 3(4vr) 3 7!' 3(4vr) 3 7! ' 3 (4tt) 3 7!' V ; 

and for a massles Dirac fermion: 

191 896 n 32 , n 40 

A = 7 , Bi = — — -r— - , B 2 = - — r^-r , and B 3 = -- — — - - . (4.12) 

3-7!' 3(4vr) 3 7!' (4tt) 3 7!' 3 (4tt) 3 7! v ; 

Thus the universal contribution (4.10) to the entanglement entropy becomes 

log(m5) for a massless conformal scalar 



S 7\ ' , . .. (4.13) 
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\og(md) for a massless Dirac fermion . 



To extract the corresponding expression for the fermion from our results, we consider 
eq. (B.41) for the universal contribution on a hyperbolic waveguide in D = 6. Replacing 
Vol(H 4 ) with Vol(§ 4 ) = 87T 3 -RVl5 and 1/R 2 -> -1/R 2 in the subsequent factors, we 
have the analogous contribution for a spherical waveguide. Then setting m = and 
a = 1, we recover precisely the universal contribution for a massless fermion given above 
in eq. (4.13). The corresponding expression for the conformal scalar can be obtained 
directly from eq. (2.49) after setting m = 0, £ = 1/5 and a = 1. Unfortunately the 
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result in this case is: S univ = — j^qq log(m<5). Hence we have a discrepancy of a factor 
of -5 between our calculation and the corresponding result in eq. (4.13)! Unfortunately, 
at this point, we have not been able to track down the source of this discrepancy despite 
checking both approaches in many ways, e.g., verifying the central charges in eq. (4.11) 
with the heat kernel techniques used here. 
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A Separation of variables in the case of spin-| fields 

The main goal of this appendix is to prove the separation of variables (2.11) in the case 
of spin-| fields. Here we focus on the waveguide geometry Ai a = C a x E> d discussed 
throughout the main text. However, in Appendix B, we also consider the hyperbolic 
waveguide C a x M d and hence we note that the present discussion is equally well ap- 
plicable to the latter case since the essential ingredient is the product structure of the 
waveguide geometry. We begin by reviewing our spinor notation. The spinors are 
associated with the orthonormal frames on A4 a 

< < 9v» = , (A.l) 

whereas the Clifford algebra in the orthonormal frame is generated by D matrices 7 a , 
satisfying the anticommutation relations 

{ 7 a , 7 b } = 2 5 ab . 

The dimension of these matrices is 2^tJ and the associated D(D 
SO(D) rotations are given by 

= 1Y1 .7 J • 

The latter satisfy the standard SO(D) commutation rules 

r aft ^cdl rfec ad s:acbd rbd ac , xad be f\ a \ 

[a , cr J = o a —da —da + o a , (A. 4) 



(A.2) 

— l)/2 generators of 
(A.3) 
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while the commutator of a with 7 C is 

[a a \ 7 C ] = 5 V - 5 ac 7 6 . (A.5) 
The covariant derivative of a spinor may be written in terms of as follows: 

V a = e£V M , V M = d„ + -o- 6c o; m6c , u^ bc = e u b (d^e cp - T^e ca ) , (A.6) 

where are the usual Christoffel symbols. With this definition, the following anti- 
commutation relations can be shown to hold [47] 

[V ai V b }ij = ~R abcd a cd ^ . (A.7) 

In particular, since the 7 a matrices are covariantly constant, one can use eqs. (A. 2) and 
(A.7) to verify the following identity for the iterated Dirac operator: 

f . f = -f 2 = -( 7 a V a ) 2 = -5 ab V a V b + I . (A.8) 

It is the heat kernel of this operator on M. a that we use in eq. (2.8) to evaluate the 
partition function in the case of spin-| fields. 

Since the waveguide geometry has a product structure, it is natural to choose the 
first two 7-matrices a = 1,2 to construct the Dirac operator on C a and then the rest 
are used to build the restriction of Dirac operator to E> d . Now, we observe that due to 
the various relations in eqs. (A. 2), (A. 4) and (A.5), the restriction of Dirac operator to 
C a anticommutes with the restriction of Dirac operator to S d . Therefore we can split 
the iterated Dirac operator into separate wave operators on C a and § d . 

It is convenient to make this discussion more explicit by choosing the following 
representation of the Dirac matrices on M. a : 

7 1 = a 1 (8) ^ , 7 2 = a 2 ® I d , (A.9) 
7 a = a 3 <g> 7 a " 2 for a = 3, • • • , D . 

where a 1 are the standard Pauli matrices while 7 a and Id is the Dirac matrices and 
unit matrix of the Clifford algebra on S d . Hence the latter matrices have dimension 
2 L rf / 2 J x 2l d / 2 l. With this explicit representation, the Dirac operator takes the form 

f Ma =fc a ®I d + a 3 ®f sd (A.10) 

where yc a an d are respectively the Dirac operators that would be constructed on 
C a and E> d alone. Similarly, the iterated Dirac operator appearing in the heat kernel 
becomes 

- flt a = -tl a ® h + h ® (-fl*) ■ (A.ll) 
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Hence it is clear that the heat kernel of the iterated massless Dirac operator on the 
waveguide C a x § d can be written as a product of the individual heat kernels, as in 
eq. (2.11). Hence the desired separation of variables is proved. 

Further we note that the spinors on _M a have dimension 2' D / 2 I = 2 2 L rf / 2 J . With 
the above construction, we see that the full trace over these spinor indices is properly 
accounted for by calculating the two spinor heat kernels on C a and S d separately and 
then simply taking their product in eq. (2.12). That is, we can treat the spinor trace 
on the right-hand side of eq. (2.12) as though we separately tracing over the spinor 
spaces of the two component manifolds. 

While the desired result has been established, let us make a few more comments. 
Recall that the construction of the spinor heat kernels E> d given in eqs. (2.52) and (2.53) 
involve a spinor matrix U which parallel propagates a spinor between the two points 
in the heat kernel. According to [32], this 'propagator' U(x,y) on § d can be formally 
written as follows 

U(x, y) = P exp i y u^ bc {t) a bc v»(t)dt , (A. 12) 

where a ab are the SO(d) generators (A. 3) for the d- dimensional Clifford algebra, now 
constructed with the Dirac matrices 7°. Above, the integration is along the shortest 
geodesic connecting x, y G S d , P is the path-ordering operator and v^(t) = dx^/dt is 
the tangent vector to the geodesic x M (t) G S d . With the representation in eq. (A. 9), 
this propagator becomes U(x,y) = I 2 ® U(x,y). Hence, given eq. (A. 11), this again 
demonstrates that the restriction of the iterated Dirac operator to C a commutes with 
the heat kernels on S d given in eqs. (2.52) and (2.53). This is, of course, necessary for 
eq. (2.11) to hold. 

Finally, let us comment that demonstrating the separation of variables in eq. (2.11) 
here crucially depends on the fact that we are studying the iterated Dirac operator 
(2.51) rather than working with the Dirac operator directly. In contrast, the eigen- 
spinors of the V on C a x S d (or H d ) cannot be obtained easily from the lower dimen- 
sional eigenspinors on the component manifolds C a and S d . While one might construct 
the heat kernel on M. a using eigenspinors, clearly such a construction would obscure 
the desired separation of variables and our discussion above demonstrates that in fact 
this approach is not needed. 

B Waveguides with hyperbolic cross-section 

Throughout the main text, we considered the waveguide geometry with a spherical 
cross-section. In this appendix, E> d is replaced by a <i-dimensional hyperbolic space M. d 
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and we show that Renyi entropy is essentially the same with an obvious sign change 
R 2 — > — R 2 in various curvature contributions. 

In the following, we mimic our previous discussion of the case of a spherical cross- 
section. In particular, massive scalars and fermions will be considered separately, and 
within each of these, odd and even dimensions d require separate consideration. Fur- 
ther, both the simple regularization of the heat kernel (2.8) and the ^-function approach 
(2.17) will be used to evaluate the Renyi entropy. Recall that while the regularization 
in eq. (2.8) reproduces the whole structure of the Renyi entropy, the ^-function method 
(2.17) eliminates the power law divergences and retains only finite and logarithmically 
divergent terms. The general remarks and formulae of section 2 do not require any 
modifications apart from obvious replacements S d — > M. d and R 2 — > —R 2 and will not 
be repeated in the following discussion. Instead we focus on the evaluation the heat 
kernel K^d and (^}. 

B.l Renyi entropy for a massive scalar 

There is a vast literature which considers the scalar heat kernel on the hyperbolic space 
M d , e.g., see [48] 

,(o) 1 ( -1 d\ n 

^ + Ahx,y) (4?rt)1/2 y 27rR 2 sinhpdf)/ 



KLAt^y) = \^^r pTp ) e-^-^, (B.l) 

K§U a (t, X ,y) = e- 3 ^* ( - ~ l - iXf^t), (B.2) 



, m2n+2 [t,x,y)-e { 2 nR 2 smhpdp) ^ 

where n is or a positive integer; p is the geodesic distance between x and y measured 
in units of R; and 

(„) -J2R f°° pe-4£ 

The scalar ^-function was computed in [31] where it was shown that for odd d > 3 

^ = R ^)Z]tm d) if^t* + 1/2,2 " * " 1/2) 1 (B ' 4) 

where B denotes the usual beta function with B(x, y) = T(x) T(y)/T(x + y) and 

b 2 = R 2 m 2 -£d(d-l) + ^Y^~ ■ (B ' 5) 

Note that b 2 here in the case of M. d is the same from its counterpart (2.26) on § d with 
the replacement R 2 — > —R 2 . The coefficients gf\ were introduced previously — see 
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eqs. (2.35) and (2.36), as well as the surrounding discussion. For even d > 4 



(B.6) 

where the coefficients h% d are given in eqs. (2.45) and (2.46) (as well as the surrounding 
discussion). In both eqs. (B.4) and (B.6), one should think of Vol(H d ) as a formal 
regulated volume for the cross-section of the waveguide geometry. While we express 
our results below in terms of this volume, it may be more appropriate to think in terms 
of the Renyi entropy density along the entangling surface. 

Odd dimensions 

Let us assume that d = 2n + 1 and take the limit of coincident points in eq. (B.l), then 

-(0), 



(d-2)/2 



2fc+l 



K^}(t, x, x) takes the following general form [49] 



K m2n+1 (t,x,x) = (47rt)n+1/2 e ^ • 



where from (B.l) it follows that P„_ 1 (x) is polynomial of degree n — 1 with rational 
coefficients. For n — 0, P_i(x) = 1 while for n > 0, P^2i(x) is polynomial of degree 
n — 1 with rational coefficients: 



n-1 
j=0 



For example, let us write out the first few polynomials 



P (0) (x) = 1 



2 



pW( x ) = 1 + - x 
3 



p(°\ x ) = l + 2x + ^x 2 , 
P 3 (0) (x) = 1 + Ax + —x 2 + —x 3 . (B.9) 

As one may surmise from these examples, a^n^i = 1 for n > 0. To simplify the 
following discussion, we also denote = for j > 0. We may evaluate the Renyi 
entropy by combining the above heat kernels (B.7) with the usual expressions given in 
the main text (i.e., eqs. (2.4), (2.8), (2.12) and (2.18)), which yields 

& = i+£ Vol(tf"«) f * fl^g . (B.10) 

24 a v J&2 t (Aixt) n+1 / 2 v ; 
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We can easily evaluate the first few divergent terms in this expression as 



- M^Volgg) 1 1 ^ 
a ' div 12 a (47r) d / 2 5 d \ d (d - 2) R 2 

Here we have used that in the vicinity of t — 0, P^l^t/R 2 ) ~ 1. As expected, we 
see that the leading divergence obeys the expected 'area law.' Continuing to higher 
orders in t, one can reconstruct all the UV divergences. We note that the structure 
of these divergences matches those for the case of E> d up to an obvious sign flip in the 
background curvature R 2 — > —R 2 , e.g., one might compare the above expression with 
n = 1 (d = 3) to eq. (2.31). 

A simple approach to evaluate the finite terms in Sa^ is to proceed in the spirit 
of dimensional regularization by setting 5 = but treating d as unspecified variable in 
the integration over t in eq. (B.10). This yields 

rf-3 

(o) _ 1 + a Vol(H d ) ^ (Q) rN _ d/2)bd -2j (BU) 

This expression is also valid for d = 1 if we omit the sum and substitute j = 0. In the 
limit mR ^> 1, the leading behaviour in eq. (B.12) is given by 

(0) 1 + a (-I)^tt / D _ 2 , D-2m D -* 

° a. finite. ~ <^a , . , D-2 ._. "^S 



ajinite ^ Q ^^p^ H 2 i? 2 

,n own J0) 



X 



(D-2)(D-3^-a^^+...j (B.13) 



where = Vol(H d ) and we have again used ciQ( D _ 3 y 2 = 1- Further we have simplified 
the above result with the following identity 

r(l - D/2)T(D/2) = — ^ = (-I)^tt (B.14) 

sin(y) 

where the last equality applies because we are only considering odd D. If we com- 
pare this result with eqs. (1.3) and (1.4), we see that the expression gives precisely the 
expected area term. We can also compare this result with eq. (2.37) for the spheri- 
cal waveguide with odd D. In this case, the second contribution above matches the 
corresponding term in eq. (2.37) with R 2 — > — R 2 if 



<4°U = ^ ~ j , (b.is) 

x > 2 1Z 
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which matches the linear coefficients in the examples given in eq. (B.9). Of course, the 
R 2 — > — R 2 behaviour is precisely that expected of a curvature contribution, as given 
in eq. (2.38). 

Alternatively we can apply the ^-function method (2.17) to derive S^ inite . Ac- 
cording to eq. (B.4), Cni (0) = an< ^ therefore, from eqs. (2.4) and (2.18), we find 

(B.16) 

.2=0 

Comparing this result to eq. (2.33), we see that the differences between the M. d and 
S d cases are not accounted for by a change in the sign of the curvature alone. While 
the first term in eq. (2.33) matches that the (entire) result above with the replacement 
R 2 —R 2 , there is also an extra integral contributing on the right hand side of 
eq. (2.33). Recall that this integral vanishes in the limit mR — > oo, however, for 
large but finite mR, it contributes exponentially suppressed terms to S^\ inite , e.g., see 
(2.32). In the main text, we speculated that these exponential contributions probe the 
topology of the waveguide geometry. In this case, it would be natural that they vanish 
here where the M. d cross-section is topologically trivial. Lastly, we note that we may 
evaluate eq. (B.16) using eq. (B.6) and compare the result with (B.13). Agreement in 
this comparison again yields eq. (B.15). 

Even dimensions 

Below we evaluate the structure of all of UV divergences for the particular case of 
d = 2. Then using ^-function approach, we evaluate (only) the universal contributions 
for any even d > 4. 



S. 



(o) 

a, finite 



1 + a 
24a 



dC {0) 
Vol(M d ) ^ 



dz 



d = 2 (D = 4): 

In this case, we substitute n = in eq. (B.2)and take the limit of coincident points 



t 



< 2 ) (t,x,x) = e-^/(° 2 )(0,t) = :: ^^ / —, - N dx (B.17) 



2AR^t Jo sinh^) 



t 

f oo 



d X e-* 2 (l-4,x 2 + ...)=—(l-^ + ... 



2ttH Jo V 6i? 2 J Ant V 12R 2 
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This result is then used as usual to evaluate the Renyi entropy, yielding 

& = ^« f f e^(v— + (B , 8) 



i±^Vol(H 2 > 
A8na 



W 2 + H " ^) l0g(m5) 



T3.19) 



Hence the leading divergence reveals the expected 'area law' behaviour, while with a—1, 
the second term is precisely the result given in eqs. (1.3) and (1.4) for D = 4. Further, 
this result agrees with the § 2 counterpart (2.42) with the usual R 2 —> —R 2 replacement. 

General even d > 4. 

While we could determine the entire pattern of divergences in the Renyi entropy in 
higher even dimensions by extending the method considered above, here we focus our 
attention on only the universal contributions. The latter are readily evaluated using 
the £- function approach (2.17) instead. 

Evaluating scalar ^-function (B.6) at z = and substituting the result into eqs. (2.4) 
and (2.17), we find that universal contribution to the Renyi entropy is given by 



(0) 1 + q Vol(W) (o) 
12a (A7r) d / 2 T(d/2)R d ^ M 



r roo x 2k+l f_yi\k+V 

-1 J BK 

(B.20) 



k=0 



4 / -= -dx- 

e 2nx + 1 k 



\og(m5) . 



Expanding in the limit mR » 1, yields 



(0) _ l + a (-I)^e / (D-2nD-3m-l)m^ N 

b °>™*>- 12a (4vr)^rp/2)V m 12 R 2 +---) [o ^ md )> 

(B.21) 

where At, = Vol(H d ). The above results are in full agreement with their counterparts for 
the spherical waveguide in eqs. (2.44) and (2.47) provided that we replace R 2 — >■ —R 2 . 

B.2 Renyi entropy for a massive fermion 

The spinor heat kernel on the hyperbolic space was evaluated in [32] as 

p ( — 1 d \ n ( p 



Ku&i (t, x, y) = U(x, y) cosh l - ( cosh " ) - — , (B.22) 

e2 n+iv , ,yj v ,yj 2\27iR 2 dcoshpJ \ 2) (int) 1 / 2 v J 

K^l (t, x, y) = U{x, y) cosh p - ( 2 ^ 2g J sh J " ( cosh (p, t), (B.23) 
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where x,y are two arbitrary points of the hyperbolic space; n = 0,1,2,...; p is the 
geodesic distance between x and y in units of R; the matrix U (x, y) is the parallel 
spinor propagator from x to y; and 



/*> (47rf)^cosh(p/2)/ Vcoshp-coshp P ' 1 J 

Alternatively, one can use ^-function method (2.17) to evaluate the partition func- 
tion Z a . The spinor ^-function on W d was computed in [32]. For odd d > 3 



fc=0 



where gj^ 2 '* are defined by 



(d-2)/2 (d— 1)/2 

II (^+/)= E ^ V "- (B.26) 

j=l/2 fc=0 



On the other hand, for even d > 2 



(1/2)^ ^VOI(H^) ( ^\ (1/2) ^ 

4 / t-^ fife 



„.2fc+l 



l)(x 2 + m 2 R 2 Y 



(B.27) 



where /iQ 1 ^ 2 '' = 1 and fojj.d for even d > 2 are defined by 



(d-2)/2 (d-2)/2 

II (^ 2 +/)= E ^M Vfc - (B.28) 

j=l fc=0 

Clearly the (^-function approach (2.17) is the most efficient way to determine the 
universal contributions in the Renyi entropy (2.4). However, as stressed previously, this 
method eliminates all of the power law divergences, while the sharp cut-off in eq. (2.8) 
allows us to keep track of the entire structure of the UV divergences appearing the 
Renyi entropy. Therefore, in what follows, we exploit both approaches to shed light on 
the structure of the Renyi entropy in the case of the spin-| field. 
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Odd dimensions 

It fol 
form 



It follows from eq. (B.22) that for d = 2n + 1, K^ 2) (x, x, t) takes the following general 



K^l(x,x,t) = ; A JS,J Ui, (B.29) 



'(1/2) f_ ^ _ Pn /2 \t/R 2 " 
(47rt)n+l/3 

where I n+ i is a unit matrix of dimension n + 1. This matrix is the remnant of the 
spinor propagator U(x,y) in the limit of coincident points. Pn 2 \x) is polynomial of 
degree n with rational coefficients 

n 

p£' 2 \x) = a f,n )xj , iorn>0. (B.30) 

3=0 

The first few of these polynomials are given by 

P^\x) = 1 , 
Pi 1/2 \x) = l + ±x, 

p£ /2 \x) = 1 + ^x + ^x 2 , (B.31) 

„fi/2i/ x 7 259 2 15 , 

P„ w ^(x) = 1 + -x H x 2 H ar . 

3 v ; 2 60 8 

By definition eig 1 ^ = 1 and as we confirm below, in general, 

o , m= (0 -2)( g -3) 

Eq. (B.29) leads to the following expression for the Renyi entropy, 
Q (i/2) _ 2 n (l + o) 2ra+1 f°° dt Pn /2 (t/ R 2 ) H 

The integral above is divergent in the vicinity of £ = 0. However, since Pn^ 2 \x) is 
a polynomial of degree n and the denominator contains a half-integer power of t, all 
the divergences are simply power-law (and thus non- universal). Using a^^f 1 = 1, the 
leading contributions take the form 

(1/2) _ (l + o) ^A^ / 1 D-2, , a l!^a\ 1 \ 

12a ( J D-2)(27r) i2 f i U fl - 2 D - 4 1 i? 2 J «^-4 + ' ' ' J ' 1 J 

where Ay, = Vo\{J§P~~ 2 ). Hence, the leading divergence exhibits the expected 'area law' 
behaviour. Continuing the process of expanding the integrand in eq. (B.33) near t = 0, 
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all power law divergences can be evaluated for any fixed value of D. We can compare 
the above expression to the results for a spherical waveguide for D = 3 and 5 given 
eqs. (2.56) and (2.64), respectively. Then up to expected substitution R 2 — > —R 2 in the 
curvature contributions, we see that the divergences in these two different geometries 
agree. 

To compute the finite part of the Renyi entropy, one can use 'dimensional regu- 
larization' approach introduced in the previous section. In this scheme we merely set 
5 = and keep n in eq. (B.33) unspecified. As a result, the power law divergences are 
eliminated and we find 



\ n (1/2) 

^finite- a 48(27T)V5F^ RV L 



\ a n (1/2) 



a 48(2^)™^ (mR) 2 i 



1 + " (-I)^tt ( rn D-2_(xm D-2m D ^ 

- 24a ( 27 r)^r(D/2)V2 ^ ^ 2 iP + " 

where we have simplified the final expression with eq. (B.14). Note that we are ex- 
panding the final result in the limit mR ^> 1. The leading contribution in this limit 
precisely matches the expected area term (1.3) with the pre-factor given by eq. (1.6). 

The sub-leading term above also agrees with eqs. (2.57) and (2.66) for D = 3 and 
5, respectively, if we replace R 2 —> —R 2 and substitute = and cl{^ = 1/2 

using eq. (B.31). We might remark that, in fact, the (^-function method can also be 
applied here to derive the same result. In particular, this calculation confirms the result 
a^pl 3 given in eq. (B.32). Following our discussion for the scalar fields, we can use 
this general result to identify the form of this subleading term for arbitrary D. We can 
then re-express this contribution in the covariant form: 

D — 2 f 

5 univ = 24 iDjermion j d D ~ 2 <7 Vh Tl(h) m D ~ 4 . (B.36) 

Again lZ(h) is the Ricci scalar of the metric induced on the entangling surface while 
the coefficient "iDjermion is given in eq. (1.6). This expression is only applicable for odd 
D > 5. 

Even dimensions 

The computation of the Renyi entropy for even dimensions is, of course, similar to 
that for odd dimensions. However, a systematic evaluation of all the divergences using 
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the simple regularization of the heat kernel (2.8) requires more effort for even d. The 
extra computational complications originate from the fact that heat kernel of the Dirac 
operator on H 2 cannot be expressed in terms of elementary functions. Instead we have 
the integral representation in eq. (B.24). 

To demonstrate the procedure, we start from the special case d = 2 and evaluate 
all the divergences in this case. The extension of this computation to higher even d 
presents no conceptual difficulties. Therefore rather than pursuing this approach (2.8) 
to construct a voluminous general expression which contains all of the (nonuniversal) 
power law divergences for general even d, we shift our focus to only evaluating the uni- 
versal logarithmic contributions. In this case, the (^-function approach (2.17) provides 
the most efficient method to produce a general result. 

d = 2 (D = A): 



In this case we need eq. (B.23) with n = 0. In particular then, the limit of coinci- 
dent points gives 

K^ /2 \t,x,x) = fg /2) (0,t)I 2 . (B.37) 

Again the two-by-two identity matrix I2 corresponds the coincident point limit of the 

spinor propagator U(x,y) on H 2 . Given the expression in eq. (B.24), we evaluate 

f (V ' 

'H 2 



/la (0, t) by expanding for small t 



tf/ 2) (M) = (A ,3 p/ - f°*e-* 2 / 4 coth( 



ytx 



dx 



(4n)*R\/i Jo v 2i? 

= i( 1 + 6^ + ---)- (R38) 

Here we only explicitly show the terms which contribute to the divergences in the Renyi 
entropy and the ellipsis denotes terms which only make finite contributions. Proceeding 
as usual with this result, we find 

s (i/2) = (1 + a) Vol(H 2 ) r at & _ tm , / + t 



a 48tt J 5 2 t 2 \ 6R 2 

1 + ol . ( 1 , 1 

487ra 



(-1 + (2m 2 - ^) log(m«J) + ...J (B.39) 



where At; = Vol(H 2 ). Of course, the leading divergence above corresponds to the usual 
'area law' term. For a = 1, the second term precisely matches eq. (1.3) with pre-factor 
given by eq. (1.6) with D = 4. Finally, eq. (B.39) agrees with the analogous result 
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(2.62) for a spherical waveguide after replacing R 2 — > —R 2 . 

General even d > 4 

Here we use the spinor ^-function (B.27) to compute the partition fucntion (2.17) 
and then the Renyi entropy (2.4). Evaluating eq. (B.27) at z = and substituting the 
result into eq. (2.17) leads to the following universal contribution to the Renyi entropy 



S. 



(1/2) (1 + a) Vol(H° 



a 12(27r) d / 2 T(d/2)R d 

(d-2)/2 



X 

k=0 



L Jo 



00 x 2k+1 , (-R 2 m 2 ) k+1 

— dx + - — 

e 2nx - 1 k + 1 



log(m5) . (B.40) 



Substituting d = 2, we recover the universal part appearing in eq. (B.39) above. Substi- 
tuting in other explicit values of d, one can generate universal contributions for higher 
even dimensions. Thus, for instance, we obtain with d — 4 (D — 6) 

3S.SL(d = 4) = -%^Vol(H 4 ) ( m * - 2^ + Iog(mtf) . (B.41) 



96vr 2 a \ R 2 30 i? 4 , 

We can also expand the above expression in the limit mR ^ 1 to determine a 
general expression for the leading terms: 

(1 + a) i^ 2 {D ~ 2)2{D ~ 3) ^ I ^ Xodm5) 

Sa —~ 12a {2^V{D/2) V 24 i? 2 + j log ^ ' 

(B.42) 

where As = Vol(HI D ~ 2 ). The leading term has the expected form (1.3) with the 
pre-factor given (1.6). The next-to-leading term reveals a new universal curvature 
contribution (1.5). Let us turn to the entanglement entropy by setting a = 1 and then 
write this subleading contributions in a covariant form as 

D — 2 f 

S univ = jDJermion / d D ~ 2 a Vh 11(h) m D ~ 4 log(m£) . (B.43) 

Again 1Z(h) is the Ricci scalar of the metric induced on the entangling surface and the 
coefficient jDjermion is given in eq. (1.6). Of course, this contribution only appears for 
even D > 4. 
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